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Abstract. We propose a many-body formalism for Cooper pairs which has similarities to the one we recently 
developed for composite boson excitons (coboson in short). Its Shiva diagram representation evidences that 
-/V Cooper pairs differ from JV single pairs through electron exchange only: no direct coupling exists due 
to the very peculiar form of the BCS potential. As a first application, we here use this formalism to derive 
Richardson's equations for the exact eigenstates of N Cooper pairs. This gives hints on why the N(N — 1) 
dependence of the A^-pair ground state energy we recently obtained by solving Richardson's equations 
analytically in the low density limit, stays valid up to the dense regime, no higher order dependence exists 
even under large overlap, a surprising result hard to accept at first. We also briefly question the BCS 
wave function ansatz compared to Richardson's exact form, in the light of our understanding of coboson 
many-body effects. 

PACS. XX.XX.XX No PACS code given 



1 Introduction 

It is commonly accepted that the Pauli exclusion principle 
plays a key role in superconductivity. None the less, the 
precise way Pauli blocking transforms a collection of single 
Cooper pairs into a superconducting condensate, still is an 
open problem. This precise understanding goes through 
the study of Cooper pairs not within the grand canon- 
ical ensemble as done in the standard theory proposed 
by Bardeen, Cooper and Schrieffer (BCS), but within the 
canonical ensemble. To handle the Pauli exclusion prin- 
ciple between a fixed number of interacting fcrmions is 
however known to be quite difficult when these fermions 
are paired. Turning to the grand canonical ensemble makes 

■ the task far easier. This is why superconductivity has been 
tackled this way, the two procedures being equally valid in 

■ the thermodynamical limit. Yet, adding fermion pairs one 
by one constitutes the one and only way to fully control 
the increasing effect of Pauli blocking from the dilute to 
the dense regime of pairs. 

Five years after the BCS milestone paper pQ on super- 
conductivity, Richardson succeeded to solve this N-body 
problem formally [2, 3 |IH[5] (see also reference [El[7]) ■ He 
showed that the exact eigenstates of the Schrodinger equa- 
tion for an arbitrary number N of pairs can be expressed 
in terms of N parameters, R±,... Rn which are solutions 
of N coupled non-linear equations, the energy of these N 
pairs reading as En = R\ + ■■■ + Rn- Although this ex- 
act form is definitely very smart, to use it in practice is 



not that easy: Except in the infinite N limit for which the 
BCS energy has been recovered [5], the equations giving 
Rn have not been, up to now, analytically solved 
for arbitrary N and arbitrary interaction strength. The 
only approaches are numerical [8ll9llT0] . This probably is 
why Richardson's solution has not had so far the atten- 
tion it deserves. Nowadays, these equations are commonly 
addressed numerically to study superconducting granules 
having small number of pairs 8 , 1T H121IT5] . 



Last year, we decided to tackle again these Richard- 
son's equations because we wanted to reveal the deep con- 
nection which has to exist between two well-known prob- 
lems, namely the one-pair problem solved by Cooper and 
the many-pair problem considered by Bardeen, Cooper 
and Schrieffer. These two problems have intrinsic simi- 
larities: In both cases, there is a "frozen" core of non- 
interacting electrons. Above this core, there is a potential 
layer with attraction between up and down spin electrons 
having opposite momenta. In the one-pair problem, this 
layer contains one electron pair only, while in the stan- 
dard BCS configuration, the potential layer is half-filled - 
a symmetrical potential on both sides of the Fermi level 
just corresponds to fill half the layer. It is clear that, by 
adding more and more pairs to the frozen Fermi sea, we 
must go from the one-pair problem studied by Cooper to 
the dense regime studied by Bardeen, Cooper and Schri- 
effer. 
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Although, at the present time, such a continuous pair 
increase does not seem easy to experimentally achieve, this 
increase can at least be seen as a gedanken experiment to 
study the evolution of the energy spectrum when the filling 
of the potential layer is changed, in order to understand 
the exact role of the Pauli exclusion principle in supercon- 
ductivity. This procedure can also be seen as a simple but 
well-defined toy model to tackle the BEC-BCS crossover 
since, by changing the number of pairs, we change their 
overlap. An overlap change has already been considered by 
Eagles Q3] , and also by Leggett [TS] , through the change of 
the interaction strength between pairs. In their approach, 
the number of pairs is fixed, so that Pauli blocking does 
not change when the overlap changes while it increases 
when the pair number increases. Consequently, the two 
types of overlap change do not involve the same physics. 
This is why to change the overlap by changing the pair 
number at constant potential is a fully relevant problem, 
complementary to the one studied in the past by Eagles 
and by Leggett. We wish to mention that Ortiz and Dukel- 
sky have applied Richardson's approach to the BEC-BCS 
crossover problem with a different perspective, the over- 
lap being varied by changing the interaction strength|16j. 
along Eagles' and Leggett's idea. Some interesting com- 
parisons between the BCS ansatz and Richardson's solu- 
tion also follow from their work. 

Since the Richardson's procedure allows one to fix the 
pair number and vary this number at will from one to half 
filling, we seriously reconsidered solving these equations 
analytically in order to better understand how supercon- 
ductivity develops from a collection of single pairs through 
the N dependence of their ground state energy. By turn- 
ing to the dimensionless form of Richardson's equations, 
we succeeded to solve these equations analytically at low- 
est order of density in the dilute limit on the single pair 
scale [17]. Indeed, these equations do have a small dimen- 
sionless parameter which is the inverse of the number of 
pairs N c from which overlap between noninteracting sin- 
gle pairs would start. This allowed us to demonstrate that, 
for N arbitrary large but N/N c still small, the energy of 
N Cooper pairs reads, at lowest order in 1/N C , as 



sity of states po, the kinetic energy of N free pairs above 
the frozen Fermi sea, is given by 



E N = N 



*F 



N — 1 

Po 



1 - 



N — 1 



Nr 



(1) 



cf is the Fermi level energy of the frozen sea. 
e c w 2f2exp (— 2/poV) is the single pair binding energy for 
a small potential amplitude V (weak-coupling limit), po is 
the density of states within the potential layer, taken as 
constant. It linearly increases with sample size. N c = poe c 
while Nq — pqQ is the number of free pair states in the 
potential layer, Q being the layer extension. 

A N(N — 1) dependence in the energy of N pairs sug- 
gests interaction treated at lowest order in density. In spite 
of it, the above result fully agrees with the textbook en- 
ergy obtained in the dense BCS configuration. In other 
words, all higher order terms in the N dependence of the 
energy cancel exactly, even under strong overlap. Indeed, 
the first term of Eq. |T]) is the exact energy of N pairs in 
the normal state whatever this pair is. For a constant den- 



E 



{normal) 
N 



2 [e Fo + (e Fo + 1/po) +■■■ +(e Fo + (N- l)/p )} (2) 

which is exactly equal to the first term of Eq. ([1]) . 

If we now turn to the condensation energy in the BCS 
configuration, obtained, for a number of pairs correspond- 
ing to fill half the potential layer, we find, according to 
Eq.© 



E 



normal 



N 



Nq e c 1 2 -2/poV /q\ 

e n = — y = 2^0^ e < pa (3) 



This result exactly matches the energy poA 2 /2 obtained 
by Bardeen, Cooper, Schrieffer within the grand canonical 
ensemble using their wave function ansatz, since the gap A 
reads as 2w c exp (— 1/poV) where 2w c is just the potential 
layer extension SI. 

The validity of a N(N — 1) interaction term over the 
whole density range seems to indicate that either Cooper 
pairs are not involved in many-body effects higher than 
2x2, or some magic cancellation takes place even in the 
dilute limit on the single pair scale. This strongly indicates 
that some unrevealed physics must hide behind such a 
surprising A-dcpcndcncc, hard to accept at first. 

As the Pauli exclusion principle is said to play a key 
role in superconductivity while Pauli blocking between 
N paired fermions is commonly known to be difficult to 
handle properly, it can appear of interest to approach 
many-body effects with Cooper pairs through a compos- 
ite boson formalism similar to the one we have success- 
fully developed for the many-body physics of excitons in 
semiconductors [IB] ■ 

The main purpose of the present work is to settle such 
a formalism. The physics being fully determined by the 
Hamiltonian, the major difference between an exciton gas 
and a set of Cooper pairs of course lies in the potential. Ex- 
citons interact via the Coulomb potential between its car- 
riers. For Cooper pairs, we here take the usual "reduced" 
BCS potential without questioning it. While its relevance 
has been proved in many physical effects, its main ad- 
vantage surely is its simplicity which makes it "solvable" 
- even more than originally thought by Bardeen, Cooper 
and Schrieffer, as seen from Richardson's works and the 
quite recent analytical solution we found to his equations. 

Coulomb potential is long-range, the BCS potential, 
taken as separable, is short range. This can be a reason 
for Cooper pairs to stay bound at large densities while ex- 
citons break through a Mott transition when the density 
increases. However, to our opinion, the crucial difference 
between excitons and Cooper pairs lies in the fact that 
usual excitons are made of fermion pairs having two de- 
grees of freedom while the electrons which interact by the 
BCS potential have opposite momenta, so that electron 
pairs have one degree of freedom only. Actually, there 
also are excitons made of pairs with one degree of free- 
dom: Frenkel excitons. Those exist in organic materials 
while Wannier excitons are found in inorganic materials. 
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The latter are made from a free electron and a free hole, 
attracted by intraband Coulomb processes. By contrast, 
Frenkel excitons are made of atomic excitations on ion 
sites which are delocalized into exciton by Coulomb pro- 
cesses between atomic levels. It is worth noting that we 
have found the same TV (TV — 1) dependence for the hamil- 
tonian mean value taken between TV ground state Frenkel 
excitons [19 , while Wannier excitons have been shown to 
have terms in N(N - 1)(TV - 2) and higher [20). 

Another important difference between excitons and 
Cooper pairs is that Cooper pairs are said to stay bound in 
the dense regime, i.e., under strong overlap, while excitons 
dissociate through a Mott transition. As a result, excitons 
when they exist, always are in the dilute limit while the 
relevant regime for superconductivity is the dense regime. 
Due to this, creation operators for single exciton eigen- 
states are relevant operators to tackle the exciton many- 
body physics while the single Cooper pair operator is 
probably not a relevant operator in the dense regime. We 
can nevertheless develop a composite boson many-body 
formalism not for correlated pairs but for the free pairs 
out of which the BCS condensate is made. This is what 
we here do. 

As a first interesting outcome of this formalism, we 
clearly see that, due to the very peculiar form of the re- 
duced BCS potential, two pairs of free electrons with op- 
posite spins and opposite momenta have an interaction 
scattering which is a succession of a fermion exchange be- 
tween pairs followed by a fermion interaction inside one 
pair. Since fermion exchange physically comes from the 
Pauli exclusion principle, this formalism evidences that 
two electron pairs interact, within the BCS potential, due 
to Pauli blocking only. However, as the Pauli exclusion 
principle acts between any number of pairs, exchange in- 
teraction scatterings - which originate from this Pauli ex- 
clusion - should a priori exist between more than two 
pairs. This strongly questions the N(N — 1) dependence 
of the TV-pair energy we found: Why Pauli-induced TVxTV 
exchanges do not show up through higher order terms in 
the energy? 

The coboson formalism is capable of exactly handling 
Pauli blocking between an arbitrary number of compos- 
ite bosons. Since the exact eigenstates of TV pairs have 
been shown to follow from Richardson's equations, a rel- 
evant first application of this formalism is to address to 
the eigenenergies of TV pairs in order to show how Richard- 
son's equations follow from this formalism, and to possibly 
understand the origin of the eigenenergy TV dependency 
we found. In doing so, we see that TV — 2 pairs stay un- 
changed when the Hamiltonian H acts on TV pairs. Since 
these pairs have one degree of freedom only, they cannot 
exchange their fermions in order to generate higeher or- 
der exchange Coulomb scattering as in the case of Wannier 
excitons. This can be the physical reason for the ground 
state energy of TV pairs to depend on TV as TV(TV — 1) only, 
with no higer order term, whatever TV, a result hard to 
accept, especially when pairs strongly overlap. 

The paper is organized as follow: 



In section [JJ we present the coboson formalism ap- 
propriate to many-body effects between the free electron 
pairs on which Cooper pairs are constructed. We derive 
the Pauli and interaction scatterings for these free pairs. 

In section [31 we use this formalism to rederive Richard- 
son's form of the exact eigenstates for TV = 1, 2, 3, • • • pairs 
interacting through the reduced BCS potential, in order 
to see how the solution for general TV develops. We then 
derive this general TV solution explicitly. 

In section III, we physically analyze the role of the 
Pauli exclusion principle in a collection of Cooper pairs. 
We, in particular, show that the l/(Ri — Rj) terms in 
Richardson's equations readily follow from Pauli scatter- 
ings for fermion exchanges between electron pairs. The 
Richardson's parameters Ri do have TV different values 
just because of Pauli blocking between Cooper pair com- 
ponents. As a direct consequence, the TV-pair ground state 
must be fundamentally different from the BCS ansatz. We 
then question this ansatz in the light of our general under- 
standing of the many-body physics of composite bosons. 
In this section, we also briefly discuss the major difference 
between Cooper pairs and Wannier excitons which could 
possibly explain why the energy of TV Wannier excitons 
has terms in TV (TV — 1)(TV — 2) and higher while they do 
not exist for Cooper pairs. 

In the last section, we conclude. 



2 Commutation technique for free fermion 
pairs making Cooper pairs 

In our recent works on the many-body physics of com- 
posite bosons - essentially concentrated on semiconductor 
excitons - we have proposed a "commutation technique" 
which allows an exact treatment of Pauli blocking be- 
tween the fermionic components of these composite bosons 
(cobosons in short). They appear through dimensionless 
"Pauli scatterings" which describe fermion exchanges in 
the absence of fermion interaction. These dimensionless 
scatterings, when mixed with energy-like scatterings com- 
ing from interactions between the coboson fermionic com- 
ponents, allow us to deal with fermion exchanges between 
any number of composite particles in an exact way. For a 
review on this formalism and its applications to the many- 
body physics of semiconductor excitons, see Refs. [18] . 

We here construct a similar formalism for the free elec- 
tron pairs on which Cooper pairs are made. 

2.1 Exchange scattering 

We consider free fermion pairs with zero total momentum 

pi = «W- k (4) 

In the case of Cooper pairs, creates a up-spin electron 

with momentum k while &l k creates a down-spin electron 

with momentum — k. These 0^ pairs have one degree of 
freedom only, namely k. This has to be contrasted to the 
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most general fermion pairs , such as Wannier ex- 

citon pairs, which have two. pairs actually have some 
similarity to Frenkel exciton pairs |19j . the index k being 
then replaced by the excited ion site n. 

It is straightforward to show that the creation opera- 
tors of these free fermion pairs commute 



f3l,Pt] =0 



(5) 



These free pairs thus are boson-like particles. It however 

fx 2 

is worth noting that while (a k ) = simply follows from 
the anticommutation of operators, the cancellation of 
(/?£) does not follow from Eq.©, but from the fact that 

4. 2 4. 2 4. 2 

(Pjj contains (a' k ) . The (p k ) cancellation which comes 
from Pauli blocking, may appear to be lost when working 
with pair operators instead of single fermion operators. 
We will see that this Pauli blocking is yet preserved in the 
commutation algebra for free fermion pairs we develop. 

If we now turn to creation and annihilation operators, 
their commutator reads 



= <5k'k - D] 



k'k 



(0) 



where the "deviation-from-boson operator" of two zero- 
momentum free fermion pairs -Dk'k reduces to 



Ac 



?k'k 



4% 



(7) 



This operator which would be zero for fermion pairs taken 
as elementary bosons, allows us to generate the dimen- 
sionlcss Pauli scatterings for fermion exchanges between 
composite bosons in the absence of fermion interaction. 
Following our works on excitons|18j. these are formally 
defined through 



A 



k'k, 



EMU^ + te^K (8) 



By noting that 



iW,0t] = <W? P = [bl k b- k ,0t] (9) 
it is then easy to show that 

Ac' lkl ,/1] = 24/ klk2 4i,k 2 (10) 

This leads us to identify the Pauli scattering of two zero- 
momentum free fermion pairs appearing in Eq. (JSj) , with a 
product of Kronecker symbols 



A 



k' 2 k 5 
k^ k. 



tfkjk/kl.k^k. 



(11) 



Such a simple expression results from the fact that these 
pairs are made of two free fermions, but also from the fact 
that they have one degree of freedom only. Actually, this 
Pauli scattering is just the one we expect for fermion ex- 
changes between (ki,k2) pairs in the absence of fermion 
interaction, as visualized by the Shiva diagram of Fig. (la). 
Indeed, from this diagram, it is clear that we must have 
(k^ = k 1; k 2 = k 2 ) and (— k 2 = — k x , — k^ = — k 2 ) which 
just gives the delta factors of Eq. flTTT) . 



-k' 



- k 'i 

k'. 



k' 



■k' 









\ / 
> 




(a) 




• 



-k 



-k 



(b) 



p 2 



-p, 



^ - k 

• — i— 



p 2 
-p. 



■p, 

p, 



+ (p\ 



(c) 

Fig. 1: Shiva diagram of free pairs 



(a) Pauli scattering A (^J for electron exchange 

between two free pairs (ki, k2), as given by Eq. (|lip . Up 
spin electrons are represented by solid lines, down spin 
elec trons by dashed lines. 

(b) | The BCS potential given in Eq. (fT4")) transforms a k 
pair into a k' pair, with a constant scattering —V, in the 
case of a separable potential «k'k = —Vw^'W^. 

(c) Interaction scattering X ( p? ) between two free 

pairs, as given in Eg. (ITS)) . Since the BCS potential acts 
within one pair only, scattering between two pairs can 
only come from exchange induced by the Pauli exclusion 
principle. 



2.2 Interaction scattering 

We now turn to the interaction scatterings resulting from 
fcrmion-fermion interaction. For a free fermion hamilto- 
nian 



ffo = £k ( a k a k + b lh 
Eq.© readily gives 

[H , Pi] =2e p /?p 



(12) 



(13) 



In standard BCS superconductivity, these fermion pairs 
interact through the reduced potential 



Vbcs = ^>k'k/3 k A 



(14) 



We will show below that this potential must be taken as 
separable «k'k = —Vw^w^ with moreover w k = w k in 
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order to possibly find the TV-pair eigenstates of Hq + Vbcs 
analytically. 

It is of importance to note that this potential funda- 
mentally is a (lxl) potential in the fermion pair subspace 
since fermion k interacts with one fermion only of the 
other species, namely fermion (— k)(see Fig. (lb)). As a 
crucial consequence, this prevents direct interaction be- 
tween two zero-momentum pairs. The only way these pairs 
feel each other, i.e., interact in the most general sense, is 
through the Pauli exclusion principle. 

For this (lxl) potential, we do have 



[Vbcs, Pi] 



(15) 



where 7 P 
potential' 



= J^k Pw v ^p while Vl, that we will call "creation 
of the free fermion pair p, is given by 



7 P K a P 



& V-p 



(16) 



The general property of creation potentials is that they 
give zero when acting on vacuum. As now shown, this 
operator allows us to generate the interactions of the p 
pair with the rest of the system. 

While the 7 P part of Eq. (fT5|) commutes with this 
is not so for the creation potential VI. Using Eq. ©, its 
commutator precisely reads 



k 



(17) 



This allows us to identify the interaction scattering for 
zero-momentum free pairs, formally defined as |18| 



X 



P2 P2 
Pi Pi 



^Pl P 2 



(18) 



with a sequence of one (2x2) fermion exchange between 
two pairs and one (lxl) fermion interaction inside one 
pair. Indeed, this sequence leads to 



P 2 P2 ^ _ 

pi pi 



E{v ik A( p k ^)+(pWpi)} 



(19) 



W P 2 .P2^Pi,Pl) ^P2,Pl 



When inserted into Eq. (TT8"]) . this readily gives Eq. (ITTl) . 

This interaction scattering is visualized by the diagram 
of Fig.(lc): the free pairs pi and P2 first exchange an 
electron. As for any exchange, this brings a minus sign. In 
a second step, the electrons of one of the two pairs, p^ or 
p' 2 , interact via the BCS potential. It is clear that, since 
the BCS potential has a (lxl) structure within the pair 
subspace, the scattering between two pairs can only result, 
as ahead said, from electron exchange between pairs, i.e., 
Pauli blocking. This diagram evidences it. 

It is worth noting that electron exchange and electron 
interaction do not play a symmetrical role in this interac- 
tion scattering. Indeed, process in which the interaction 



takes place before the exchange 
Fig.lc - would lead to 



instead of after as in 



£ A (pl" 



«k P] 



— ^p' 1 ,P 2 ^P 2 ,P2' !; P2Pl 



(20) 



which is definitely different from the first term of Eq.(19). 

In the next section, we use this commutation formal- 
ism to rederive the equations that Richardson has ob- 
tained for the eigenstates of N Cooper pairs through a 
totally different route. The new derivation we have pro- 
posed, through its diagrammatic support, enlightens some 
important physical aspects of this exactly solvable prob- 
lem. 



3 Richardson's equations for N Cooper pairs 

In order to better grasp how these equations develop, we 
are going to increase the number of pairs in the potential 
layer one by one, starting from a single pair. 



3.1 One pair 

We first consider a state in which one free pair (k, — k) 
is added to a "frozen" Fermi sea \Fq), i-e., a sea which 
does not feel the BCS potential. This means that the iVk 
prefactors in Eg. (TUT) cancel for all k belonging to \Fq) in 
order to have Vbcs \Eq) = 0. 

Note that this "one-pair" state actually contains iVo + 1 
electron pairs, Nq being the number of pairs in the frozen 
sea, so that this state is a many-body state already, but 
in the most simple sense since the Fermi sea \Fq) is just 
there to block states by the Pauli exclusion principle. This 
Fermi sea mainly brings a finite density of state for all 
states above it, a crucial point to have a bound state in 
3D whatever the weakness of the attracting BCS potential. 

By choosing the zero energy such that Hq \Fq) = 0, 
Eqs. (|13ll5l) gives the hamiltonian H = H + Vbcs acting 
on a one-free-pair state as 



HBt\F ) = [H,Bt 
We then note that 



IF, 



o) = (2e k 8t + 1 t + v£)\F ) (21) 



v£ \F ) = 



(22) 



since the v\^ p factor included in the 7 k part of V k , brings 
v *P a i a P \ F o) = Ukp&lp&_p \F ) = 0; 

Next, we subtract E\0^ \Fq) to the two sides of the 
above equation, with E\ yet undefined, but assumed to 
be different from any 2et- We then divide the resulting 
equation by (2ek — E\). This gives 



(H - E 1 ) 



1 



2e k - E_ 



-Pi Wo) = Pi \F ) + 



1 



2e k - E 1 



7 k 1^0} 

To go further and possibly obtain the one-pair eigen- 
state of the hamiltonian H in a compact analytical form, 
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it is necessary to approximate the BCS potential by a sep- 
arable potential u kp = —Vw k w p . The operator 7^ in Eq. 
(fT5l) then reduces to 



where (ft is given by 



(24) 



(25) 



The interaction part of H acting on two free pairs is 
visualized by the diagram of Fig. [2j In Vkik 2 )> one P a i r 
stays unchanged while the other pair suffers a BCS inter- 
action. In |W / kik 2 ); the two pairs exchange an electron and 
then one pair interacts. These diagrams evidence the fact 
that, due to the (lxl) structure of the BCS potential, two 
pairs can interact by fermion exchange only as a result of 
the Pauli exclusion principle. 



If we now multiply Eq. (1231) by Wk and sum over k, we end 
with 



{H - E 1 )B\E l )\F Q ) 



2e k - E 1 



ft \F ) 



where the operator B'(E) is defined as 



(26) 



Pi (27) 



Eq. (j2"6")) readily shows that B^ (E\) \Fq) is one-pair eigen- 
state of the hamiltonian H with energy E\ , provided that 
the bracket in the RHS is zero, i.e., E\ fulfills 



2e k - E x 



(28) 



This is just the well-known equation for the single pair 
energy derived by Cooper. 



3.2 Two pairs 



We now add two pairs to the frozen sea \Fq). Eas. (|13ll5p 
yield 

HPlA. \F ) = ([H,Pl] Pl+Pl [H,0l]) \F ) 

= (2^ + 2^)4X1^) (29) 
+ |^ kl k 2 ) + |W klk2 ) 

The last two terms come from interactions between the 
four electrons of the two pairs. The first term of Eq. fj 15[) 
readily gives |Vk 1 k 2 ) as 

|Vi Clk2 )=( 7kl 4 2 +7iX)l F o) 

= -^(^k 1 /3 k2 +</3 kl )/3 t |^o) 
The second term of Eq. (TIB"]) yields 

\W klk2 }=(vl0l+ptvl)\F o ) 

To calculate it, we again use commutators. Since acting 
on the frozen sea \Fq) gives zero (see Eq. (|2"2"]l ). we find 
from Eqs. fiMty 



(30) 



(31) 



i^o = E^S)ttifo) 



(32) 



" k 'i -§-- 



k\ 



k. 



+ (k,- 



Fig. 2: Shiva diagram for the interaction part |Vk lk2 ) + 
Wkik 2 ) of the Hamiltonian H acting on two free pairs, as 
given in Eqs.(30)and (32) 

To go further, we subtract £3 \F ) to the two 
sides of Ea. ([2T))) . with E 2 yet undefined. We split E 2 as 
R\ + R 2 and we multiply the resulting equation by 
w kl w k2 / (2e kl - i?i) (2e k2 - R 2 ). This yields 



{H - E 2 )Bt x {R x )Bt^R 2 )\F Q ) = 



^k 2 /3 k2 



Vwl 



2e k — R- 



+ (1«2) \F ) 



2V 



ki 



fit 

2 (2e kl -J?i)(2e kl - #2) kl 



ft\F ) (33) 



the last term coming from the exchange interaction term 

I WW 

As a last step, we sum over (ki , ka). The sum over (lt 2 ) 
in the first bracket readily gives 



2e k - i?2 



ft 



(34) 



To calculate the sum over (ki,k2) in the second bracket, 
we first note that 



1 



(2e kl - Ri) (2e kl - R 2 ) 
1 1 



(2e kl - Ri) (2e kl - R 2 ) 



1 



(35) 



(Ri - R 2 ) 



which is valid, provided that R\ ^ R 2 , a condition that 
we can always enforce since the unique requirement is to 
have i?i + R 2 = E 2 . For u£ = u> k , we then find 



kik 2 



(2e kl -i?i)(2e kl - R 2 )^ 
1 



2V5 klkl w kl /3t i fi \F ) 



(Ri - R2) 



B<<(Ri)-Bi(R 2 )} (36) 
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Summation over (ki,k2) of Eq. (1331) then yields 



(H - E 2 )B^{R 1 )B^(R 2 )\F Q ) = 



1 - V 



2V 



2e k — i?2 R\ ~ R2 



+ (1 O 2 
^\F ) (37) 



The above equation evidences that B< \R\)B' \R%) \Fq) 
is two-pair eigenstate of the hamiltonian H with energy 
E2 = R\ + R2 provided that the bracket in the above 
equations is zero, i.e., (Ri, R2) fulfill two equations, known 
as Richardson's equations for two pairs 



2V 



2e k — i?i R\ — R-2 



= (l*+2) 



(38) 



3.3 Three pairs 

We now turn to three pairs in order to see how these 
equations develop for an increasing number of pairs. Two 
usually is not generic, while three most often is. We will 
here see that when H acts on three pairs, one at least 
among the three pairs stays unchanged. This is a step 
toward understanding why the N dependence of the TV- 
pair ground state energy is in N{N — 1) only, with no 
term in N(N — 1)(N — 2) and higher, as the validity of 
our low density result extrapolated to the high density 
BCS regime, seems to indicate. 
We start with 



HPtAjl \ F o) 

+pIpI [h^I]}\f ) 



{ 



(39) 



Using Eqs. (|13ll5[) . we again split the above equation into 
a kinetic part and two interaction parts 

H PlAA a l F °> = + 2£ k, + 2e k3 ) PlAA 3 l F °> 
+ |^k lk2 k 3 ) + |Wk lk2 k 3 ) 

(40) 

As in the case of two pairs, the BCS potential generates 
direct processes which are given by 

iv klk2k3 > = (iIpIpI+iIpIpI +7^4 A) i F °> 

(41) 

since 7^ and ftl commute. It also generates exchange pro- 
cesses, which appear as 

iw klk2k3 > = (vIpIAs +pI v IpI 1^0) 

To calculate them, we again use commutators and the fact 
that V£\F ) =0. Eq. (JS) allows us to rewrite |W klk2 k 3 ) 



in a more symmetrical form as 

{ [Vl > Pl 2 ] Pi + Pi [Vl , Pt 3 ] + Pi [Vl , Pi] } \F ) 



(43) 



The interaction part of Eq. (|4*Tj| . namely |Vk 1 k 2 k 3 ) + 
|Wkik 2 k 3 ), is represented by the diagrams of Fig|3] 
|Vkik 2 k 3 ) bas interactions inside a single pair, two pairs 
staying unchanged. |W / kik 2 k 3 ) contains processes in which 
the pair suffering the potential exchanges one of its elec- 
trons with a second pair, the third pair staying unchanged. 
There are three similar contributions, obtained by circular 
permutations. 



k' • 



" k 3 

K 
-k 



k 



A 



k\ 



k 3 
" k 3 

K 

-k 



" k , 

k. 



Fig. 3: Shiva diagram for the interaction part |Vkik 2 k 3 ) 
tU klk2k3 ) of the Hamiltonian H acting on three pairs 



Using Eq. ([T9"]) for the interaction scattering, we end 
with 



|^ kl k 2k3 ) = 2V(6 klk2 w kl (3l(3l + 2 permit \F ) (44) 

which has close similarity with |W kl k 2 ) given in Eq. ((32|) . 
To go further, we insert Eq.([ST]) for U klk2k3 ) and Eq.([ii|) 

for |W kl k 2 k 3 ) into (gDJ); we subtract E 3 plAA 3 \F Q ) to 
both sides, with E3 written as R\ + R2 + R3, and we mul- 
tiply the resulting equation by 

w kl Wk 2 w k3 / (2e kl - i?i) (2e k2 - R 2 ) (2e k3 - R 3 ). This 
yields 



(H - £ 3 )<(i?i)<(i? 2 )<(i? 3 ) \F ) 
{Bl(Ri)Bl(R 2 ) 
+2 perm.} |F ) 



-2V 



'k lk2 ^k! 



{< {R3) [(2e kl -R 1 )(2e kl -R2)^ 
+2 perm.}/? 1, \F ] 



(45) 



8 



Monique Combescot, Guojun Zhu: Coboson formalism for Cooper pairs used to derive Richardson's equations 



We then sum over (ki, k2, Its). By calculating the sums of 
the two brackets as for two pairs, Eqs. (l34|) and (|36f then 
yield 

(H - E 3 )B\R 1 )B^(R 2 )B\R 3 ) l^o) - 
{B\R 2 )B\R 3 ) 



2V 



+ 



2V 



2e kl — ^1 Ri — Ri R3 ~ Ri 

+ 2 perm.}/3 f |F ) (46) 



This leads us to conclude that the three-pair state 
B^(Ri)B^(R 2 )B^(R 3 )\F ) is eigenstatc of the hamilto- 
nian H with the energy E3 = R\ + R2 + R3, provided 
that (i?i, i?2, R3) fulfill the three equations, 



wl 


2V 


2V 


2e k - Ri 


R\ — i?2 


R\ — JI3 


wl 


2V 


2V 


2e k - R2 


h — — - 

R2 — R3 


i?2 — Rl 


wl 


2V 


2V 


2e k - Rs 


R 3 — Ri 


R3 — R2 



(47) 



3.4 N pairs 



The above commutation technique can be easily extended 
to N pairs. As nicely visualized by the diagrams of Figs|2] 
and [3J the effect of the BCS potential on these N pairs 
splits into direct and exchange processes: In the direct set, 
one pair is affected by the (lxl) scattering while the other 
N — l pairs stay unchanged. In the exchange set, this pair, 
before interaction, also exchanges one of its electrons with 
another pair, the remaining N— 2 pairs staying unchanged. 
This understanding shows that an increase of pair num- 
ber above two, does not really change the structure of 
the equations since N — 2 pairs stay unchanged, the pair 
which exchanges its fermions with the pair suffering the 
interaction being just one among (N — 1) pairs. 

The procedure is rather straightforward once we have 
understood that either (N — 1) or (N — 2) pairs stay un- 
affected in the BCS interaction process. The general form 
of the iV-pair eigenstates ultimately appears as 



(h-e n )b^(r 1 )---bHr n )\f ) = o 



(48) 



with En = Ri + h -Rjvj these Rn's being solutions of 

N coupled equations 



i = ^E 



E 



2V 



2e k — Ri tt? Ri ~ Rj 



for i = (1,...,N) 



Let us explicitly derive these N equations following the 
procedure we have used for three pairs. 



\Fo) 



(50) 



The hamiltonian acting on N pairs can be written in 
terms of commutators with any of these N pairs as 

= ~H,Pl]0l 2 ---Pt N \Fa) 
+ --- + 4 1 ---/3 ki _ i [H,^ 

We then use Eqs. (fnTT5|) to replace H, /3 k . by 2e ki /3 k . + 
7 k . + Vjl. . The first two contributions commute with the 

other /3 k 's, so that, when inserted into the above equation, 
they yield 



N 



2(e kl 



- e k iV )n/ 3 kj^o)+i^ kl ... kJV ) 



(51) 



where the direct interaction part is given by 



N 

2—1 m^i 

which is similar to Eq. (|4"Tj). 

The part with the creation potential Vl, is more cum- 
bersome. We again calculate it through commutators. Let 
us consider one term. We start as 

viPi i+ ,---Pi N \F ) 

= Vl^t i+1 ]pt i+2 --Pt N \Fo) 

+ ■■■ + ■ ■ ■ Wl > &] \ p o) 



+ 



(53) 



Vj^ , /? k . makes appear the exchange interaction scatter- 
ing between the k^ and kj pairs, so that the above term 
generate the exchange interaction scatterings between the 
ki pair and all the kj pairs with i < j < N. When inserted 
into Eq. (f5U|) . we find all possible exchange interactions be- 
tween fermion pairs, so that the set of terms with creation 
potential ultimately gives 

i^ k ,.. kN ) = E^AEx(g k k :) n H f o) 



P1P2 



i<3 



(54) 

which is similar to Eq. (|4"3"f . If we now use Eq. (|52"]l for the 
sum over (p' 1; p' 2 ), we end with 



\W kl ... kN ) = 2V^[5^ 2 w*X 2 '■'&„+'■ 
which also reads in a compact form as 



l*b> (55) 



(49) |W kl ... kw )=2V r /3t 



I^O) 

(56) 
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All this leads for En written as R\ 



(H-E N )(3l...(3l N \F ) 

= [(2e kl -R 1 ) + ... + (2e k2 - R N )]pl ■ ■ ■ (3^ \F ) 



(57) 



To go further, we do as before: we multiply both sides of 
the equation by w Kl ■ ■ ■ w kjv /(2e kl - Ri) ■ ■ ■ (2e k]v - Rn) 
and we sum over (ki, ■ • • , kjv). The LHS readily gives 



(H-E N )B\R 1 )---B^(R N )\F ) 



The first term in the RHS, which comes from the free pair 
kinetic energy, yields 

[(£ wkA)*'^) • • • fit (^) + • • ■] l^o) 

N 



=j9 t EII Bt ( i ^)i i? o> 

The first interaction term, induced by direct processes 
within one pair, readily leads to 



A' 



^Ec-^E- 



. 2£k s - R, 



R N to 4 Physical understanding 

4.1 Richardson equations and the Pauli exclusion 
principle 

The above derivation of Richardson's equations makes crys- 
tal clear the parts in these equations which are directly 
linked to the Pauli exclusion principle between fermion 
pairs through electron exchanges. 

From a mathematical point of view, the link is rather 
obvious: In the absence of terms in V/(Ri — Rj), the N 
equations for Ri would reduce to the same equation (|28[) . 

so that the solution would be R^ — E\ for all i. The fact 
that the energy of N pairs differs from N times the single 
pair energy E\ entirely comes from the set of (Ri — Rj)'s 
different from zero. 

Physically, the fact that En differs from NE\ comes 
from interactions between Cooper pairs. Due to the (lxl) 
structure of the BCS potential within the pair subspace, 
interaction between pairs can only be mediated by fermion 
exchanges as evidenced from the diagram of Fig. (ITc|) . 
Consequently, interactions between pairs are solely the re- 
sult of the Pauli exclusion principle. This Pauli blocking 
mathematically appears through the various <5 P ' P factors 

(5g) in the Pauli scatterings A( p ? ^ 2 ). These S factors are the 

ones of the |Wki— kjv) term. They ultimately lead to the 
various (Ri—Rj) differences in the Richardson's equations, 
as easy to follow from our procedure. 

In short, the Kronecker symbols in the Pauli scatter- 
ings of fermion pairs take care of states which are ex- 
cluded by the Pauli exclusion principle. They induce the 
2V/(Ri — Rj) terms of the Richardson's equations which 
) [[ B\R m ) \F ) (60) ultimately makes the energy of N pairs different from the 

energy of N independent pairs. 



(58) 



while contributions coming from exchange interaction pro- 
cesses appear as 



4.2 Excitons versus Cooper pairs 



2Vft 



E Sk ^ 



k, 



Lkik 2 



flt 

(2e kl -i?i)(2e k2 -i? 2 )^ kl 
B^(R 3 )---B\R N )\F ) 



An important feature of the ground state energy En for N 
pairs that this derivation possibly explains, is the fact that 
(61) the part of En coming from interaction, namely En—NE\ 
depends on N as N(N — 1) only, with no higher order 
dependence. Indeed, Eq.(l) also reads 



By using Eq. (1361) for the above bracket, we end with 
(H-E n )B^R 1 )---B^Rn)\F ) = 



N 



2V 



Ri — Ri 



2—1 k j^i 

Y,BHRrn)\F ) (62) 

This evidences that B^(Rx) ■ ■ ■ B^(R N ) \F ) is N -pair eigen- 

state of the hamiltonian H, with energy En = Ri H h 

Rn provided that all the brackets in the above equation 
cancel. These just are the N Richardson's equations writ- 
ten in Eq.(49). 



E 



N 



NEi + N{N - 1) — + 



1 



Pa 



- 



(63) 



In order to have terms in the energy in N(N — l)(Af — 
2), we need topologically connected diagrams between 3 
pairs. The diagram of Figj3] shows that, when H acts on 
three pairs, one out of them do not participate to the 
scattering, so that the three pairs are not connected. In 
the case of pairs, this is (N — 2) out of the N pairs 
which are not connected. Connections thus seem to exist 
between two pairs only. 

This however is not enough to explain that higher or- 
der terms do not exist in the energy because the energy 
of N Wannier excitons has terms in N(N — 1)(A — 2) and 
higher [5D] while the hamiltonian acting on N Wannier ex- 
citons also leaves N — 2 excitons unchanged. Indeed, let 
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B\ be the creation operator of the i exciton with energy 
E u i.e., (H - Ei)B\ |0) = 0. We do have, as in Eq. O, 



hb] b) B\ 10) = 

l\ 12 *3 I / 



{ 



H.Bl 



B f B 1 + B ] 



12 »3 
t 



if, BJ 



5;' 



(64) 



H,B 



}|0) 



To calculate it, we introduce the i exciton creation oper- 
ator V} defined as 



H,B\ 



EiB] + V} 



(65) 



which has similarity with Eqs. (|13I15I) . This operator is 
such that V} |0) = 0, as readily seen from the above equa- 
tion acting on vacuum. From it, we construct the interac- 
tion scatterings of two excitons through 



y^B]\=^{i^BiBi 



(66) 

which is similar to Eo. lfTg)) . When used into Ea. ljM)) . this 
yields 

(H-E n -E V2 -E l3 )BlBlBl\0) 

=$>tx{£(:i?K (e?) 



As in Eq. (|43|) , one out of the three excitons seems to stay 
outside in the scattering process. 

Here comes the crucial difference between Wannier ex- 
citons and Cooper pairs. Wannier exciton, made of 
pairs, have two degrees of freedom. As a direct conse- 
quence, two electrons and two holes can be associated in 
two different ways to form two excitons. It is possible to 
show [THI that 



B]B\ 



VA(" j )BlBl 



(68) 



where A ( " 3 ) is the Pauli scattering of two excitons, de- 
fined, as for /?£ pairs, through 



B m , B\ 



D R' 



(69) 



The fact that, in the first term of the RHS of Eq.(67), 
the B\ exciton does not participate to the scattering, is in 
fact somewhat artificial because, using Eq. (|68|) . we could 
as well write this first term as 



E 



BtBtBl{X(ll)+X{ll)}^ m l) (70) 



This shows that the ?2j *s) excitons can actually be 
involved in 3 x 3 connected diagram, as shown in Fig. 01 



m 



Fig. 4: Shiva diagram for energy-like exchange interaction 
between three excitons 



Let us return to the pairs making Cooper pairs. We 
see that, since the k electron with up spin is associated to 
the — k electron with down spin only, an equation similar 
to Eq. fTjSl) does not exist. As a result the k3 pair in Eq. (|4"3"| 
cannot be mixed with the kJ, pair as in Eq. (|70l) to generate 
3x3 connected diagram, as the one of Fig. |4l 



4.3 Richardson's exact eigenstate versus BCS ansatz 

The Richardson's procedure we have here rederived, gives 
the exact form of the Hq + Vbcs eigenstates as 



Bt(ij 1 )...flt(fl JV )| Fo ) 



(71) 



with B^R) given by Eq. (p?T|) . The fact that, by construc- 
tion, all the i?i's are different in order for the l/(f?; — 
Rj) factors in the Richardson's equations not to diverge, 
strongly questions the standard BCS ansatz for the N- 

pair wave function since this ansatz reduces to (f?^) |io) 
when projected into the TV-pair subspace. In this ansatz, 
all the pairs are taken as condensed into the same state. 
This is physically hard to accept for composite bosons due 
to Pauli blocking between pairs which makes each added 
pair necessarily different from the previous ones, due to 
the fact that more and more states are occupied already. 

There were in past several discussions about differ- 
ences and similarities between BCS ansatz and Richard- 
son's exact solution, or more generally a Bethe ansatz like 
Hi B[ |0) from various perspectives and physical situation 
21,22 . However, the discussions essentially focus on re- 
covering the correct energy or some other physical quanti- 
ties, not the wave function itself, more difficult to experi- 
mentally evidence. This wave function actually is attached 
to the picture people commonly have of superconductiv- 
ity. This is why a correct wave function is importance for 
physical understanding, at least. 

To discuss this problem, let us again start with two 
pairs. In a previous work[23j, we have shown, that the 
Richardson's parameters for two pairs read as R\ = R + 
iR' and R2 = R — iR' with R and R' real. In the large sam- 
ple limit, i.e. for 1/po small, the dominant terms of R and 
R 1 are given by R ~ e c + l/ po + e c / PqS7 and R' w ^2e c /po. 
By writing B^Rx) as [B^(R)+B^(Rx)-B^(R)] and sim- 
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ilarly for ST(i? 2 ), we get, from eq (12"T)) . 

S t (i? 1 )St(ii 2 ) - [B\R)] 2 = R' 2 [C|C1 - 2£t (R)D^\ 

(72) 

where we have set 



^ t = E 



^ 

(2e k -ii) (2e k -i?±iR') k 
iu k 



(2e k -i?) (2e k - i?) 2 + i?' 2 



4 



(73) 
(74) 



Eg. (|T2"|) shows that, in order to possibly replace 
B^{R 1 )B^{R 2 ) by [B f (i?)] 2 as in the BCS ansatz, we 
must neglect terms in R' 2 , i.e., in 1/po- However, these 
1/po terms are precisely those which make E\ different 
from E 2 /2 « _Ei + 1/po + e c j p^fl; so that the replace- 
ment of B> (Ri)B^ (R 2 ) by a "condensed two-pair state" 

(Bi(E 2 /2)) with E 2 different from E x is fully inconsis- 
tent because, in this two-pair operator, we would keep 
contributions in 1/po which are as large as the ones we 
drop by neglecting the RHS of Eq. (f72")) : two pairs do not 
condense into the same state. 

Actually, it is claimed that the BCS ansatz is valid in 
the thermodynamical limit when N is very large. It is pos- 
sible to show that, for N large but still in the dilute regime 
on the single pair scale, the Rfs stay two by two complex 
conjugate, the imaginary part of Rfs getting larger and 
larger as */ Ne c / po when N increases. By using a similar 
procedure as the one we used for N = 2, we hardly see 
how, starting from the exact form of the iV-pair eigenstate 
B^(Ri) ■ ■ ■ B^(Rjy) \F Q ), we can possibly recover the BCS 
ansatz with the same creation operator for all pairs when 
N ventures outside the dilute limit because nothing spe- 
cial happens in the behavior of the Rfs when N crosses 
N c . 

In a recent work, Ortiz and Dukclsky[l6] have also 
considered Richardson's equations in the thermodynami- 
cal limit. While they do recover the energy obtained from 
the ansatz, they conclude, like us, that Richardson's ex- 
act wave function is substantially different from the BCS 
ansatz in many ways. 

We wish to stress that, to the best of our knowledge, 
derivations of the validity of the BCS ansatz for the ground 
state of N pairs mainly concentrate on the energy it pro- 
vides (see, e.g., [53] and references therein). We of course 
agree that the BCS ansatz gives the correct ground state 
energy for N pairs because the energy obtained using 
this ansatz is just the one we derived from the exact 
Richardson's procedure, extrapolated outside the dilute 
limit. However, agreement on the energy by no mean proves 
agreement on the wave function. Many examples have 
been given in the past with wave functions very different 
from the exact one, although giving correct energy. Di- 
rect experiments supporting the form of the ground state 
wave function however seems to be even harder to achieve 
than the ones possibly checking the N dependence of the 
ground state energy given in Eq.(l). Nevertheless, it seems 
to us highly desirable to carefully reconsider "agreement 
with experiments" in the light of the exact Richardson's 



wave function. It is still a rather intriguing question to un- 
derstand why the minimization of the hamiltonian mean 
value calculated with this ansatz, leads to exactly the 
same energy as the one we derived by analytically solv- 
ing Richardson's equations in the dilute limit. 

It is worth noting that the reduced potential used in 
standard BCS superconductivity has the great advantage 
to allow an analytical resolution of the iV-body Schrodinger 
equation - which is quite infrequent. It however is clear 
that this potential is highly simplified. A certain amount 
of corrections are necessary to make this potential more re- 
alistic. These are going to destroy the possibility to get the 
eigenstates analytically. However, since the BCS ansatz for 
the wave function with all the pairs condensed into the 
same state - which is commonly considered as one of the 
essential features of superconductivity - has been worked 
out within this reduced BCS potential, a precise compar- 
ison between this conventional ansatz and the exact solu- 
tion of the model in the canonical ensemble, is definitely 
quite relevant to better understand the deep physics hid- 
den in this ansatz. 

Finally, we wish to stress that the possible replacement 
ofSt(i? 1 )..._Bt( J R JV ) |F ) by (B*) N \F ) is crucial to sup- 
port the overall picture of superconductivity commonly in 
mind, with all the pairs in the same state, "as an army of 
little soldiers, all walking similarly" . This picture actually 
seems a rather naive extrapolation to composite bosons, 
of the standard Bosc-Einstcin condensation demonstrated 
in the case of elementary bosons. It is hard for us to ac- 
cept that, in the case of composite bosons, Pauli blocking 
between fermionic components is not going to destroy nice 
harmony in this "army" . More work on the validity of the 
BCS ansatz in the thermodynamical limit in the context 
of the coboson nature of Cooper pairs, seems a necessity 
to more deeply understand some unrevealed aspects of ba- 
sic superconductivity as the ones at the origin of Eq.(l) 
for the iV-pair ground state energy. The coboson many- 
body formalism we have here constructed, should appear 
as quite valuable because it gives a fresh view to this fa- 
mous field, its Shiva diagram representation helping to 
support physical understanding. 



5 Conclusion 

We have constructed a coboson formalism for the electron 
pairs on which Cooper pairs are made. It has similarity 
with the one we have constructed for composite boson 
excitons. This formalism evidences that the scatterings of 
two zero-momentum electron pairs in the BCS potential, 
are mediated by the Pauli exclusion principle. No direct 
process exists 

As a first application, we here rederive Richardson's 
equations for the exact eigenstates of N Cooper pairs. 
This derivation allows us to trace back the physical origin 
of the various terms. In particular, we clearly see that N 
pairs differ from N independent pairs, due to Pauli block- 
ing only. This Pauli blocking enforces the Ri parameters 
of Richardson's equations to be all different. As a direct 
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consequence, the exact wave function for N interacting 
pairs is definitely different from the BCS ansatz, although 
the N-pair energy this ansatz provides, is the correct one 
in the large sample limit. 

The diagrammatic representation of our derivation also 
shows that, because electron pairs with zero total momen- 
tum have one degree of freedom instead of two, they scat- 
ter within the BCS potential through (2 x 2) exchange in- 
teraction scatterings only. This possibly explains why the 
./V-pair ground state energy that we have recently found in 
the dilute limit, has interaction terms in N(N — 1) but not 
in N(N — 1)(N — 2) and higher, as expected for TV-body 
problems. 

One of us (M.C.) wishes to thank the University of 
Illinois at Urbana-Champaign, and Tony Leggett in par- 
ticular, for enlightening discussions during her invitation 
at the Institute for Condensed Matter Physics where most 
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the manuscript. 
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